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The efficiency at maximum power (EMP) for tight-coupling molecular motors is investigated based 
on the constitutive relation between the flux and thermodynamic force. The motors are classified into 
four generic types (linear, superlinear, sublinear, and mixed types) according to the characteristic 
of the constitutive relation. Correspondingly, the EMP (77*) for these four types of motors is proved 
to satisfy 77* = 1/2, < n* < 1/2, 1/2 < 77* < 1, or < 77* < 1, respectively. For the given reduced 
chemical potential Au which is the ratio of the released free energy of fuel in each motor step to the 
energy scale of thermal motion at the physiological temperature, the lower and upper bounds of EMP 



are found to be [1 - W(e 1_AM )]/A/j and [W(e 



1]/ Au, respectively, where W(.) represents 



the Lambert's W function. A phase diagram with phase boundary S = 2/ Au — l/(e AA1,/2 — 1) is 
constructed, which shows how Au and the load distribution factor 5 influence on the EMP: The 
value of EMP is smaller than 1/2 when the parameter pair (Au, 8) takes value in the region above 
the boundary while larger than 1/2 below the boundary. This fact implies that the motors using 
ATP as fuel under the physiological condition can work at maximum power with higher efficiency 
if S < 0.1, which provides a thermodynamic interpretation to several experimental observations on 
kinesin and myosin motors. 

PACS numbers: 05.70.Ln, 05.40.-a 



Introduction. Molecular motors are special proteins 
which take charge of the long-distance transport in liv- 
ing cells 1]. The energetics of motors has attracted much 
attention from physicists since 1990s 041] ■ It is highly 
convinced that the operations of molecular motors abide 
by some economic principles of energetic managements 
because the motors are products of the long-term evolu- 
tion and natural selection. 

One of the possible economic principles might be that 
molecular motors work at the region of maximumpower. 
To clarify this principle, Schmiedl and Seifert Q opti- 
mized the power of molecular motors based on a generic 
model and found that the efficiency at maximum power 
(EMP) for tight-coupling molecular motors is 1/2 when 
the motors work in the linear non-equilibrium region 
where the released free energy of fuel in each motor step 
is very small. They also deduced that the EMP can be 
larger or smaller than 1/2 when the motors operate in 
the region far from equilibrium [9J. These results have 
been well confirmed and developed by other model sys- 
tems [IM1, which can be regarded as the counterparts 
of the EMP for heat engines investigated by many re- 
searchers [I6I - I25I ]. 

Recently, Van den Broeck et al. derived the upper 
and lower bounds of EMP for tight-coupling molecular 
motors, and found that these two bounds can be reached 
when the transition state is close to the reactant state or 
the product state, respectively. However they expressed 
the bounds of EMP as the functions of thermodynamic 
force which depends on both the released free energy of 
fuel in each motor step and the external load on the mo- 
tors. Both of them are in fact not independent of each 
other when the motors work at maximum power. Can 



we express the bounds as the explicit functions of the 
released free energy of fuel in each motor step? In addi- 
tion, most of the above researches reveal that the EMP of 
tight-coupling molecular motors can be larger or smaller 
than 1/2, which depends on the load distribution fac- 
tor (i.e., the position of the transition state) and the 
released free energy per fuel molecule (for example the 
hydrolysis energy of ATP). How do these quantities in- 
fluence on the EMP of molecular motors? In this work, 
we will address these problems by investigating the EMP 
for tight-coupling molecular motors. The motors are clas- 
sified into four generic types (linear, superlinear, sublin- 
ear, and mixed types) according to the characteristic of 
the constitutive relation between the flux and thermo- 
dynamic force. We obtain the corresponding ranges of 
EMP for these four types of molecular motors and ex- 
press the bounds of EMP as the explicit functions of the 
released free energy of fuel in each motor step. A phase 
diagram is construct, which implies that the motors us- 
ing ATP as fuel under the physiological condition can 
work at maximum power with higher efficiency for the 
small load distribution factor, which provides a thermo- 
dynamic interpretation to several experimental observa- 
tions on kinesin and myosin motors. 

Model and optimization. We consider a simple transla- 
tional molecular motor with equivalent discrete sites X n 
(n = 0, ±1, ±2, • • • ) with distance I between the nearest 
neighbor sites proposed by Schmiedl and Seifert As 
schematically shown in Fig. [T] the motor will consume 
one fuel molecule and move against the external load / 
in each step. That is, the chemical step and the mechani- 
cal step are tightly coupled. This process can be regarded 
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FIG. 1. (Color online) Scheme of the energy landscape V(x). 
I and 8 represent the step size and the position of transition 
state, respectively. 



as a chemical reaction 



Fuel + X n =± X n+ i + Products, 



(1) 



where 



= fc+e s f l and cj_ = fe_e^ are the for- 



ward and backward rate constants, respectively. In the 
expressions of rate constants, we have not explicitly writ- 
ten out the energy scale (fc^T) of thermal motion at the 
physiological temperature. Thus fl in fact represents 
fl/ksT in the present work. k + (fc-) depends on the 
bare forward (backward) rate constant and the concen- 
tration of fuel (products). The load distribution factor 
< S < 1 indicates the position of the transition state, 
i.e., the distance between the transition state and the re- 
actant state. The extreme cases 5 = and 1 correspond 
to the situations that the external load merely influences 
on the backward or forward rate constants, respectively. 
Thermodynamic consistency [9j requires 



uj + /cj_ = c 



(2) 



where A/i is called the reduced chemical potential which 
represents the ratio of the released free energy of per fuel 
molecule in the above reaction to the energy scale (fc^T) 
of thermal motion at the physiological temperature. 

Let us consider the steady state where each sites have 
equal probability. The entropy production rate may be 
expressed as 



a = (lu + — cj_ ) ln(cj + /oj- ) 



(3) 



with omitting an unimportant prefactor. The net flux is 
defined as 



J 



(4) 



Thus the fundamental thermodynamic relation (entropy 
production rate = flux x thermodynamic force) enlight- 
ens us that the conjugate thermodynamic force can be 
expressed as 



F = ln(w+/u;_) =A(i-fl 



(5) 



with considering Eqs. ©-(UJ). Then the constitutive re- 
lation between the flux and the thermodynamic force can 



be expressed as 



e SF_ e -(l-S)F 



(6) 



The energy input per unit time for the tight-coupling 
motor can be expressed as G = (/cbTA/i) J. The power- 
output is the difference between the energy input and the 
energy dissipation per unit time, which can be expressed 
as 

P = G- k B Ta = k B T(An - F)J. (7) 

Thus the efficiency can be defined as 

r, = P/G = 1 - F/Ajj,. (8) 

It is stressed that we adopt the traditional definition of 
efficiency in the present work rather than the Stokes ef- 
ficiency, rectification efficienc y o r sustainable efficiency 



22|. 



proposed in recent Refs. 

Our discussions so far can essentially be regarded as 
the concrete representations of the general considerations 
for the EMP of molecular motors in Refs. [HI 



ii, m 



Now maximizing the power with respect to the external 
load /, we obtain 



J / J' + F = A/i 



(9) 



from Eq. ([7]), where J' represents the derivative of the 
flux J with respect to the thermodynamic force F. Sub- 
stituting the above equation into Eq. ([8]), we derive that 
the EMP satisfies the following relation 



1 



1 + J'F/J' 



(10) 



J 



,'superlinear linear,-' 








FIG. 2. Schematic diagram of four generic types of consti- 
tutive relations between the flux J and the thermodynamic 
force F. 

Bounds of EMP. As we have done for heat engines in 
Ref. [261 ] . we may classify motors into four generic types 
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according to the characteristic of constitutive relation be- 
tween the flux J and the thermodynamic force F. As 
shown in Fig. [2J the motor is of linear type when the 
constitutive relation is linear, i.e., J increases uniformly 
with F. Similarly, a motor is of superlinear (sublinear) 
if the increasing rate of J is enlarged (reduced) with the 
increasing of F in the whole curve. Additionally, the mo- 
tor is of mixed type if the increasing rate of J is enlarged 
with the increasing of F in some segment while reduced 
in the other one. In fact, the motor of this type crosses 
over from the sublinear type to the linear type and then 
to the superlinear type with the increasing of F, or vice 
versa. 

The linear, superlinear and sublinear types of con- 
stitutive relation can be mathematically expressed as 
J = J'F, J < J'F, and J > J'F for any finite F, re- 
spectively. Considering < F < A/i and Eq. (fTO)). we 
can easily derive 



V* = 1/2, 
< 77, < 1/2, 
1/2 < rj, < 1, 



linear tpye; 
superlinear type; 
sublinear type. 



(11) 



The constitutive relation of mixed type is a little compli- 
cated, which can be mathematically expressed J > J'F 
in some segment (sublinear type) while J < J' F in the 
other one (superlinear type), and J = J' F at the demar- 
cation point (linear type) between both segments. Cor- 
respondingly, the EMP can be larger or smaller than 1/2, 
or equal to 1/2 for the motor of mixed type, depending 
on which region F is located in. That is, the EMP for 
the motor of mixed type satisfies < 77* < 1 . 

Now we discuss three specific examples as follows. The 
first one is the situation that A/i <C 1 so that F C 1. 
In this case, J is approximately a linear function of F. 
According to Eq. (TIT)) , we derive 77* — > 1/2 for A/i -C 1, 
which confirms the previous result obtained in Refs. I9| — 

The second one is the limit case of 8 = 1. On the one 
hand, substituting 8=1 into Eq. ([6]), we derive the flux 
J = fc_(c F — 1) which displays the superlinear behavior 
(,/ < J'F) for any finite F. Thus Eq. |Qj} implies < 
77* < 1/2. On the other hand, substituting the flux into 
Eq. ([5]), we can obtain the optimized F„ = A/i — 1 + 
W(e 1_A/1 ) where W(.) is the Lambert's W function [27, 
[28j . Substituting the optimized F* into Eq. (|5J), we finally 
drive 



v*V = i) = 



1 - W(e 1 " A ^) 
A/i 



(12) 



This analytic expression is plotted in Fig. [31 which reveals 
that the EMP indeed satisfies < 77, < 1/2. 

The third one is the other limit case of 6 — 0. On the 
one hand, substituting 6 = into Eq. ((6]), we derive the 
flux J = fc_e AAl (l — e~ F ) which displays the sublinear 
behavior (J > J'F) for any finite F. Thus Eq. ([IT]) 



implies 1/2 < 77, < 1. On the other hand, substituting 
the flux into Eq. ([9]), we can obtain the optimized F* = 
A/i + 1 - W(c 1+Af *). Substituting the optimized F* into 
Eq. (|HJ), we finally drive 

W(e 1+A ^) - 1 



770 = r^(S = 0) 



A/i 



(13) 



This analytic expression is plotted in Fig.[3j which reveals 
that the EMP indeed satisfies 1/2 < 77* < 1. 




FIG. 3. Curves of EMP for tight-coupling molecular motors. 
771 and 770 represent the EMP for 5 = 1 and which can 
be analytically expressed as Eqs. (|12[1 and {T3j, respectively. 
The triangles (a), squares (■), stars (-fc) and circles (o) cor- 
respond to the EMP for 8 = 0.1, 0.3, 0.5 and 0.7, respectively. 

In general case, from Eq. ([6]) we can prove that the 
motor is of superlinear type (J < J'F) for 1/2 < 8 < 1. 
Thus the EMP should be smaller than 1/2 for 1/2 < 8 < 
1. On the other hand, it is not hard to find that the 
motor is of mixed type for < 6 < 1/2 through simple 
calculations from Eq. ([6]), and that the corresponding 
constitutive relation displays the behavior similar to the 
solid curve in Fig. [5J Then the EMP can be larger or 
smaller than 1/2 for different values of A/j. A natural 
problem is to investigate the exact bounds of EMP for 
given value of A/x. 

Substituting Eq. ([6]) into Eq. ([9|), we derive 

We cannot solve the analytic solution to this equation if 
5 or 1 . However, we readily see that 8 is a monoton- 
ically increasing function of F for given A/i because we 
find 88/ dF > from the above equation. Equivalently 
speaking, F is also the monotonically increasing function 
of (5 for given A/j. Additionally, Eq. ([8]) implies that the 
efficiency is the monotonically decreasing function of F. 
Thus 77* should be the monotonically decreasing function 
of 8 for given A/j,. It follows that 77* is bounded by rji 
and 770, that is 

(15) 



1 - Wie 1 -^) W(e 1+A ' i ) 
- < V* < - 



A/i 



A/i 
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which is the first main result in this work. We note that 
Van den Broeck et al. also derived the bounds of 77* 
in Ref. [Hj]. However they expressed the bounds as the 
functions of F which depends on both A/i and /. Here 
we express the bounds as the explicit functions of single 
variable A/i. 

We also find numerical solutions to Eq. (fl4l) for 8 = 0.1, 
0.3, 0.5 and 0.7 with various values of A/i. Then the 
numerical relations between 77* and A/i are achieved from 
Eq. (jHJ and shown in Fig. [3] as triangles, squares, stars 
and circles for 6 = 0.1, 0.3, 0.5 and 0.7, respectively. 
Obviously, all numerical data are indeed located in the 
region bounded by 771 and 770, which is consistent with 
inequality ()15[) . In addition, 77* is always smaller than 
1/2 for 8 = 0.5 and 0.7 while can be larger or smaller 
than 1/2 for 8 = 0.1 and 0.3. This fact is consistent with 
our above discussions because the constitutive relations 
corresponding to 8 = 0.5 and 0.7 are of superlinear type 
while those corresponding to 8 — 0.1 and 0.3 are of mixed 
type. 

Phase diagram. There are two parameters in the 
present model: one is the load distribution factor 8; an- 
other is the reduced chemical potential A/i. We will con- 
cern how these two parameters determine 77* to take value 
larger or smaller than 1/2. Since we have proved 77* = 1/2 
in the trivial case of A/i = 0, we only discuss the case of 
A/i > in the following contents. To address this prob- 
lem, we first find the condition (the relation between 8 
and A/i) to make 77* = 1/2. It follows that J'F/J = 1 
and F = A/i/2 from Eqs. © and (fT0]l, Considering 
Eq. ©, we derive 



2 

A/i 



1 



3A/V2 _ 1 ' 



(16) 



That is, when 8 and A/i satisfy the above relation, 7/* 
should always be 1/2. Since we have proved that the 
efficiency is a monotonically decreasing function of 8 for 
given A/i in the above discussion, we can readily deduce 
77* > 1/2 if 8 < 2/A/i - l/(e A ^ 2 - 1) and vice versa. 
This is the second main result in the present work. 

We draw the phase diagram in Fig. [4] according to 
the above analysis. The phase boundary is described by 
Eq. (fl6|) . When 8 and A/i take values in the region above 
the boundary, the motor can work at maximum power 
with lower efficiency (< 1/2). On the contrary, when 8 
and A/i take values in the region below the boundary, 
the motor can work at maximum power with higher effi- 
ciency (> 1/2). In Table HI we list the numerical results 
of 77* for several values of 8 and A/i and the positions 
of the pair (A/i, 8) in the phase diagram. These numeri- 
cal results indeed support our above predictions from the 
phase diagram. 

In Ref. [9( , Schmiedl and Seifert argued that the long- 
term evolution and natural selection might shape the 
ability that motors can work at maximum power with 
higher efficiency. Their insight can be easily understood 





1 1 1 
TJ<T]<M2 


1/2<7 7 <^ ^ 





A/., 



FIG. 4. (Color online) Phase diagram. The solid curve is 
the phase boundary described by Eq. (|16p which separates 
the parameter plane into two regions. The value of EMP is 
smaller than 1/2 when the parameter pair (A/i, 8) takes value 
in the region above the boundary while larger than 1/2 below 
the boundary. 



with the aid of the phase diagram according to which 
the parameter pair (A/i, 8) should be located in the re- 
gion below the boundary. For the motors using ATP 
as fuel, A/i m 20 under the physiological condition, we 
calculate 8 < 2/A/i - l/(e A ^/ 2 - 1) « 0.1 which is in 
good agreement with several experiments on kinesin and 
myosin motors where the small 8 (< 0.1) was observed 
for the main motor step 29l-32j|. 



TABLE I. The numerical results of 77* for some values of 8 



and A/i. 


A/i 


8 


Position 


V* 


5.0 


0.311 


On the boundary 


0.50 


10.0 


0.193 


On the boundary 


0.50 


20.0 


0.100 


On the boundary 


0.50 


5.0 


0.7 


Above the boundary 


0.27(< 1/2) 


10.0 


0.5 


Above the boundary 


0.20(< 1/2) 


15.0 


0.2 


Above the boundary 


0.33(< 1/2) 


5.0 


0.1 


Below the boundary 


0.65(> 1/2) 


15.0 


0.05 


Below the boundary 


0.78(> 1/2) 


20.0 


0.01 


Below the boundary 


0.85(> 1/2) 



Conclusion and discussion. In the above discussions, 
we investigate the EMP for tight-coupling motors. It 
is found that this efficiency is bounded by 771 = [1 — 
W(e 1 " A ^)]/A/i and 770 = [W(e 1+A ^) - l]/A/i as shown 
in Fig. [3l We obtain a phase diagram (Fig. [4j of param- 
eters 8 and A/i and determine the phase boundary by 
Eq. (|16|) . Considering the argument that motors might 
operate at maximum power with higher efficiency [H , we 
derive 8 < 0.1 for the motors using ATP as fuel un- 
der the physiological condition, which provides a ther- 
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modynamic interpretation to several experimental obser- 
vations on kinesin and myosin motors. Finally, we clarify 
several points as follows that have not been emphasized 
in the above discussions. 

i) It is a little subtle to take 1/2 rather than other val- 
ues as the threshold to distinguish the higher efficiency 
from the lower one for motors when we construct the 
phase diagram. The most reasonable philosophy comes 
from the Darwinism. If the efficiency of some motors is 
smaller than 1/2, then the energy dissipation is larger 
than the useful energy That is, most of the input energy 
is wasted by the motors. This kind of motors might be 
unfavorable for the natural selection, therefore they will 
easily be extinct in the evolution. Combining this consid- 
eration, the phase diagram in this work, and the exper- 
imental data on kinesin and myosin motors, we propose 
to refine and develop Schmiedl and Seifert's argument 
into an economic principle of energetic managements: the 
long-term evolution and natural selection shape the good 
performance of molecular motors such that they can op- 
erate at maximum power with efficiency larger than 1/2. 

ii) The analysis in the present work might be ex- 
tended to discuss the bounds of EMP for information 
machines or the Feynman ratchet investigated in recent 
literature (33H35|. At a glance, the EMP of the Feyn- 
man ratchet derived in our previous work [35j j might in 
fact be the lower bound because it corresponds to the 
case of 5 = 1. However, in the recent work [36j, Van den 
Broeck and Lindenberg found that the EMP for classical 
particle transport and the EMP of the Feynman ratchet 
share the same expression, which seems to imply that the 
EMP of the Feynman ratchet should be independent of 
5. It is indeed not hard to verified this point based on 
the discussions in the present work. 

iii) All discussions in the present work are merely fo- 
cused on a minimal model with discrete sites where the 
internal states are neglected. While the main prediction 
6 < 0.1 of this minimal model is in good ag reement with 
most of experimental observations Mill, a transition 
state in the range S ~ 0.3 — 0.65 has also been extracted 
[37j j based on a seven-state model. It is necessary to in- 
vestigate how the model with internal states influences 
on the bounds and phase diagram in this work. 
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